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Topological quantum states with non-Abelian Fibonacci anyonic excitations are widely sought after for the 
exotic fundamental physics they would exhibit, and for universal quantum computing applications. The frac¬ 
tional quantum Fiall (FQFi) state at filling factor v = 12/5 is a promising candidate, however, its precise nature 
is still under debate and no consensus has been achieved so far. Fiere, we investigate the nature of the FQH 

V = 13/5 state and its particle-hole conjugate state at 12/5 with the Coulomb interaction, and address the issue 
of possible competing states. Based on a large-scale density-matrix renormalization group (DMRG) calcula¬ 
tion in spherical geometry, we present evidence that the essential physics of the Coulomb ground state (GS) at 

V = 13/5 and 12/5 is captured by the fc = 3 parafermion Read-Rezayi state (RRs), including a robust exci¬ 
tation gap and the topological fingerprint from entanglement spectrum and topological entanglement entropy. 
Furthermore, by considering the infinite-cylinder geometry (topologically equivalent to torus geometry), we ex¬ 
pose the non-Abelian GS sector corresponding to a Fibonacci anyonic quasiparticle, which serves as a signature 
of the RRs state at 13/5 and 12/5 filling numbers. 


Introduction .— While fundamental particles in nature are 
either bosons or fermions, the emergent excitations in two- 
dimensional strongly-correlated systems may obey fractional 
or anyonic statistics [1, 2]. After two decades of study [3- 
13], current interest in exotic excitations focuses on states 
of matter with non-Abelian quasiparticle excitations [14—16], 
and their potential applications to the rapidly evolving field of 
quantum computation and cryptography [17-22]. So far the 
most promising platform for realization of non-Abelian statis¬ 
tics is the fractional quantum Hall (FQH) effect in the first 
excited Landau level, and two of the most interesting exam¬ 
ples are at filling factors ^ = 5/2 and 12/5. The v = hj2 
state is widely considered to be the candidate for the Moore- 
Read state hosting non-Abelian Majorana quasiparticles [14- 
16]. Experiments have revealed that the 12/5 state appears 
to behave differently from the conventional FQH effect [5, 8 ], 
and may also be a candidate state for hosting non-Abelian ex¬ 
citations. However, the exact nature of the FQH 12/5 state is 
still undetermined due to the existence of other possible com¬ 
peting candidate states. 

Several ground-state (GS) wavefunctions have been pro¬ 
posed [16, 25, 29-32] as models for the observed FQH ef¬ 
fect at ^ = 12/5 [5, 8 , 13]. The most exciting candi¬ 
date is the k = i parafermion state proposed by Read and 
Rezayi (RR3)[16]. This RR 3 state describes a condensate of 
three-electron clusters that forms an incompressible state at 
V = 13/5 [16]. One can also construct the particle-hole part¬ 
ner of the RR 3 state to describe the 12/5 FQH effect. Be¬ 
sides the RR 3 state, some competing candidates for v = 13/5 
or 12/5 exist: a hierarchy state [26, 27], a Jain composite- 
fermion (CF) state [28], a generalization of the non-Abelian 
Pfaffian state by Bonderson and Slingerland (BS) [29, 30], 
and a bipartite CF state [31, 32]. So far, the true nature of 
the 12/5 and 13/5 FQH states remains undetermined. The 
main challenges in settling this issue are the limited compu¬ 
tational ability and the lack of an efficient diagnostic method. 
For example, from exact diagonalization (ED) calculations in 
the limited feasible range of system sizes, it is found that the 


overlaps between the Coulomb GS at z/ = 12/5 and different 
model wavefunctions are all relatively large [16, 30], while 
the extrapolated GS energies of the RR 3 and BS states are 
very close in the thermodynamic limit [30, 33]. Taken as a 
whole, previous studies have left the nature of the Coulomb 
GS at zz = 13/5 and 12/5 unsettled. 

Recently, there has been growing interest in connecting 
quantum entanglement [34-37] with emergent topological or¬ 
der [38, 39] in strongly interacting systems, which offers a 
new route to identification of the precise topological order 
of a many-body state. Although characterization of entan¬ 
glement has been successfully used to identify various well- 
known types of topological order [40-47], application of the 
method to a system with competing phases still faces chal¬ 
lenges when ED studies suffer from strong finite size effects, 
and other methods such as quantum Monte-Carlo suffer from 
sign problems. The recent development of the high efficiency 
density-matrix renormalization group (DMRG) in momentum 
space [44, 54] allows the study of such systems in sphere and 
cylinder geometries, both of which can be used to make con¬ 
crete predictions of the physics of real systems in the ther¬ 
modynamic limit. Here we combine these advances, and use 
these two geometries to address the long-standing issues of 
the EQH at zz = 12/5 and 13/5. 

In this paper, we study the EQH at zz = 12/5 and 13/5 
filling by using the state-of-the-art density-matrix renormal¬ 
ization group (DMRG) numerical simulations. By studying 
large systems up to N^. = 36 on spherical geometry, we es¬ 
tablish that the Coulomb GS at zz = 13/5 is an incompress¬ 
ible EQH state, protected by a robust neutral excitation gap 
A„ Ri 0.012(e^/ZB). Crucially, we show that the entangle¬ 
ment spectrum (ES) fits the corresponding SU (2)3 conformal 
field theory (CET) which describes the edge structure of the 
parafermion RR 3 state. The topological entanglement entropy 
(TEE) is also consistent with the predicted value for the RR 3 
state, indicating the emergence of Eibonacci anyonic quasipar¬ 
ticles. Moreover, we also perform a finite-size scaling analysis 
of the GS energies for zz = 12/5 states at different shifts cor- 
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responding to the particle-hole-conjugate of the RR 3 state, the 
Jain state and BS state. Finite-size scaling confirms that the 
ground state with topological shift S = —2(3) (where RR 3 
state is expected to occur) is energetically favored in the ther¬ 
modynamic limit. Finally, to explicitly demonstrate the topo¬ 
logical degeneracy, we obtain two topological distinct GS sec¬ 
tors on the infinite cylinder using infinite-size DMRG. While 
one sector is the identity sector matching to the GS from the 
sphere, the new sector is identified as the non-Abelian sector 
with a Fibonacci anyonic quasiparticle through its character¬ 
istic ES and TEE. Thus we establish that the essence of the 
EQH state at v = 13/5 is fully captured by the non-Abelian 
parafermion RR 3 state (and by its particle-hole conjugate at iz 
= 12/5) and show that it is stable against perturbations as we 
change the Haldane pseudopotentials and the layer width of 
the system. 

Model and Method. — We use the Haldane representation 
[26, 48, 49] in which the We electrons are confined on the sur¬ 
face of a sphere surrounding a magnetic monopole of strength 
Q. In this case, the orbitals of the n-th LL are represented as 
orbitals with azimuthal angular momentum —L, —E-l-1,..., L, 
with L = Q + n being the total angular momentum. The total 
magnetic flux through the spherical surface is quantized to be 
an integer A/, = 2L. Assuming that electron spins are fully- 
polarized and neglecting Landau-level mixing, the Hamilto¬ 
nian in the spherical geometry can be written as: 

X] {rni,m2\V\m3,m4)dl^^dl^.^drn3drn^ 

mi+m2=m3+m4 

where (dm ) is the creation (annihilation) operator at the 
orbital m and V is the Coulomb interaction between electrons 
in units of /Ib with Ib being the magnetic length. The two- 
body Coulomb interaction element can be decomposed as 

2L I 

m= — l 

where {L, i; L, j\l, m) is the Clebsch-Gordan coefficients and 
V"(/) is the Haldane pseudopotential representing the pair en¬ 
ergy of two electrons with relative angular momentum 2L — I 
in n-th LL [26, 71]. Eor electrons at fractional filling factor n. 
Ns = i/~^Ne — S, where S is the curvature-induced “shift” 
on the sphere. 

Our calculation is based on the unbiased DMRG method 
[50-55], combined with ED. The (angular) momentum-space 
DMRG allows us to use the total electron number Ne and the 
total z-component of angular momentum 
good quantum numbers to reduce the Hilbert subspace dimen¬ 
sion [54]. Here, we report the result atv = 13/5(12/5) with 
electron number up to Wg = 36(22) by keeping up to 30000 
states with optimized DMRG, which allows us to obtain accu¬ 
rate results for energy and the ES on much larger system sizes 
beyond the ED limit (N^^ = 24(16) at i/ = 13/5(12/5)). 

Groundstate Energy, Energy Spectrum and Neutral Gap .— 
We first compute the GS energies for a number of systems up 
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FIG. 1: (a) The groundstate energy per electron (blue dots) corre¬ 
sponding to the V = 13/5 state. The blue line shows the extrapo¬ 
lated values obtained using a quadratic function of 1/Ne. The red 
dots shows the rescaling energy by a renormalized magnetic length 
and the red line is the linear fitting, (b) The neutral gap A„ for 
13/5 state as a function of the 1/We (6 is the layer-width parame¬ 
ter [71]). Inset: Energy spectrum versus total angular momentum 
for Ne = 21. A„ is defined as the energy difference between 
the lowest energy state (in = 0) and the first excited state (in 
/ 0 ). 


to Wg = 36 at ^ = 13/5, with a shift 5 = 3 consistent with 
the RR 3 state. As shown in the low-lying energy spectrum 
in the inset of Eig. 1(b) obtained from ED for Ng = 21, the 
GS is located in the = 0 sector and is separated from 
the higher energy continuum by a finite gap, which signals an 
incompressible EQH state. The extrapolation of the GS en¬ 
ergy to the thermodynamic limit can be carried out using a 
quadratic function of 1/Ng (blue line), or a linear fit in 1/Ng 
(red line) after renormalizing the energy by ^j2Qv/Ne to take 
into account the curvature of the sphere [56], as shown in Eig. 
1(a). We obtain the Eo/Ng = —0.38458(24) (blue line) and 
—0.38487(9) (red line), which demonstrates consistency be¬ 
tween the two extrapolating schemes. 

We also calculated the neutral excitation gap A„ at ^ = 
13/5 [57]. This is equivalent to the energy difference be¬ 
tween the GS and the “roton minimum” [5 8-60] as illustrated 
in the inset of Eig. 1(b). The roton minimum corresponds to 
the lowest excitation energy of a quasielectron-quasihole pair 
[60]. Eig. 1(b) shows A„ as a function of 1/Wg, where the 
large-system results indicate that the neutral gap approaches 
a nonzero value A„ Ri 0.012 ± 0.001 for Ng > 21. Since 
the hamiltonian in this paper is particle-hole symmetric, the 
neutral gap at 1 / = 12/5 and 13/5 are expected to be iden¬ 
tical [61]. In addition, if the effect of finite layer-width is 
considered[71], the neutral-excitation gap is reduced but still 
remains consistent with a nonzero value (Eig. 1(b)). 

Competing states. — In Eig. 2, we compare the GS ener¬ 
gies per electron of three known candidates for 1 / = 12/5: the 
particle-hole conjugate of the RR 3 state with a shift S = —2, 
the non-Abelian BS state with 5 = 2 [29], and Jain state with 
5 = 4. We find that the lowest-energy state for the Jain state 
shift (5 = 4) in larger system sizes has a total angular mo¬ 
mentum 0 , indicating that it represents excitations 

of some other incompressible state rather than the Coulomb 
GS at iz = 12/5 [32]. Secondly, the GSs with the RR 3 and 
BS shifts continue to have = 0 for the systems that 
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FIG. 2: Finite-size extrapolation of the ground-state (GS) energies 
for different shifts corresponding to different candidate states at i/ = 
12/5. All energies have been rescaled by the renormalized magnetic 
length. The angular momentum of the GS is shown whenever it is 
nonzero / 0). 


we have studied, and the extrapolation based on the result for 
10 < fVe < 22 leads to Eo/Ng = —0.3425 for the RR 3 state 
and Eq/Ns = —0.3410 for the BS state, respectively. Com¬ 
pared to the previous studies [30, 33], the extrapolation errors 
are reduced by the inclusion of larger system sizes obtained 
using DMRG. Our calculations suggest that the RR 3 state with 
shift S = —2(5 = 3) is energetically favored as the GS at 
zz = 12/5(13/5). Our results are consistent with the inter¬ 
pretation that the RR 3 state describes the true GS (see the full 
evidence below), while the other states at nearby shifts corre¬ 
spond to states with quasiparticle or quasihole excitations. 

Orbital ES. — Li and Haldane first established that the or¬ 
bital ES of the GS of FQH phase contains information about 
the counting of their edge modes [36, 39]. Thus, the orbital 
ES provides a “fingerprint” of the topological order, which 
can be used to identify the emergent topological phase in a 
microscopic Hamiltonian [36, 41-44]. 

As a model EQH state, the RR 3 parafermion state can be 
represented by its highest-density root configuration pattern 
of “ 1110011100 ... 11100111 ”, corresponding to a general¬ 
ized Pauli principle of “no more than three electrons in five 
consecutive orbitals” [62-64]. Consequently, the orbital ES 
depends on the number of electrons in the partitioned subsys¬ 
tem [71]. In Eig. 3, we show the orbital ES of three dis¬ 
tinct partitions for system size Ng = 36 for Coulomb GS. 
Eor 3n electrons in subsystem (Eig. 3(a)), the leading ES dis¬ 
plays the multiplicity-pattern 1 , 1 ,3, 6,12 in the first five an¬ 
gular momentum sectors = 0,1, 2,3,4. Eor 3n -F 1 or 
3n + 2 electrons in subsystem (Eig. 3(b-c)), the ES shows the 
multiplicity-pattern of 1, 2, 5,9 in the AL = 0,1, 2, 3 momen¬ 
tum sectors. The above characteristic multiplicity-patterns 
of the low-lying ES agree with the predicted edge excitation 
spectrum of the RR 3 state obtained either from its associated 
CET, or the “< 3 in 5” exclusion statistics rule [71, 72]. 

In addition, we vary the Haldane pseudopotentials V^(l) 
and V^(3) (keeping all others at their Coulomb-interaction 
values), and map out an ES-gap diagram which illustrates 
the robustness of the EQH state as the interaction param- 




FIG. 3: (a-c) The low-lying orbital ES of Ng = 36 are shown 
for three different partitions. The lower ES level counting in the 
sector AL^ = 0, 1,2, 3, 4 are labeled by color, where = 

— L^^^rnin with as the quantum number where the pri¬ 

mary field occurs. The entanglement gap of orbital ES of Ng = 24 is 
shown for partition (d) with 3n electrons and (e) with 3n-Fl electrons 
in the subsystem as a function of pseudopotential V^(l)/Vcoui(l) 
and V^(3)/Vco„i(3), where VcouiiO are the Coulomb values of 
pseudopotentials. The black point corresponds to the Coulomb point. 


eters are changed[65-68]. In Eig. 3, we plot the en¬ 
tanglement gap (for the lowest-L^ ES level) [36, 54] as 
a function of Vi(l)/V^o„;(l) and Vi(3)/V^o^;(3), where 
^Coii/(0 '^he Coulomb values of pseudopotentials. We 
find that the entanglement gap is robust in a region cen¬ 
tered at an approximately-fixed V^(l)/V^(3) ratio (indi¬ 
cated by the white line). Away from that, for the regime 
V'(l)/V^„,;(l) < 0.92 and Vi(3)/V^„„,(3) > 0.98, we 
find a rapid drop of the entanglement gap indicating a quan¬ 
tum phase transition. We have also studied the effect of 
the ES of modifying the Coulomb interaction with a realistic 
layer width (b) [71], and find that the RR 3 state persists until 
b/ls ^ 2 , which is qualitatively consistent with the results of 
varying V^(l) and V^(3). 

Topological Entanglement Entropy. — Eor a two- 
dimensional gapped topologically-ordered state, the 
dependence of the entanglement entropy Sa { Ia ) of the 
subsystem A on the finite boundary-cut length Ia has 
the form Sa { Ia ) = odA — 7, where TEE 7 is re¬ 
lated to the total quantum dimension 27 by 7 = In 27 
[34, 35]. We have extracted the TEE using our largest 
system, Ng = 36 [71]. The TEE obtained was 

7 = 1.491 ± 0.091, consistent with the theoretically- 
predicted value 7 = ln27 = In \/5(l + ifP') ~ 1.447 for 
the RR 3 state, where each non-Abelian Eibonacci anyon 
quasiparticle contributes an individual quantum dimension 
dp = (j) = (v/5 + l)/2 {(j) denotes the Golden Ratio). 
The appearance of dp = ^ is a signal of the emergence 
of Fibonacci anyon quasiparticles, and arises because two 
Fibonacci quasiparticles may fuse either into the identity 
or into a single Fibonacci quasiparticle [47]. This exotic 
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property makes Fibonacci quasiparticles capable of universal 
quantum computation [17]. 



FIG. 4: (a-b) The low-lying orbital ES of 'I'l and by setting 
Ly = 241s. |4'i(,^)) denotes the GS with identity 1 (Fibonacci 
0) anyonic quasiparticle, (c) Energy difference and (d) entropy dif¬ 
ference between |'l/i) and obtained from infinite DMRG on 

cylinder geometry with varying Ly. The error bars are determined 
based on results from ten different infinite DMRG calculations for 
each sector[71]. 


Topological Degeneracy on the infinite cylinder .— 
Topologically-ordered states have characteristic GS degen¬ 
eracies on compactified spaces. To access the different 
topological sectors at i/ = 13/5, we implemented the infinite- 
size DMRG in cylinder geometry with a finite circumference 
Ly [44, 69, 71]. For each value of Ly, we repeatedly calcu¬ 
lated GSs using different random initializations for the infinite 
DMRG optimization. We found that each infinite DMRG 
simulation converged to one of the two states; j'l'i) and 
These states are distinguishable by their orbital ES as shown 
in Fig. 4; |T'i) has the same ES structure as in Eig. 3(a-c), 
which matches the identity sector with root configuration 
“... 0111001110 ...”. On the other hand, |T'^) shows the 
ES multiplicity pattern 1,3,6,13,..., which identifies the 
spectrum as that of the Eibonacci non-Abelian sector with 
root configuration “... 1010110101...” [71]. Eurthermore, 
these two groundstates are indeed energetically degenerate, 
with an energy-difference per electron of less than 0.0002 
with Ly = 2AIb, while the entropy difference between 
these two states is around AS” Ri In^ Ri 0.48, consistent 
with the quantum dimension of the Eibonacci quasiparticle. 
Combining this with the fivefold center-of-mass degeneracy, 
we have obtained all the 10 predicted degenerate RR 3 GSs on 
infinite cylinder (or torus). 

Summary and discussion. — We have presented what we 
believe to be compelling evidence that the essence of the 
Coulomb-interaction ground states at 1 / = 13/5 and 12/5 is 
indeed captured by the parafermion fc = 3 Read-Rezayi state 
RR 3 , in which quasiparticles obey non-Abelian “Eibonacci- 
anyon” statistics. The neutral excitation gap is found to be a 
finite value A„ r; 0.012e^/ZB in the thermodynamic limit. 
Results for the entanglement spectrum “fingerprint” and the 
value of the topological entanglement entropy show that the 
edge structure and bulk quasiparticle statistics are consistent 


with the prediction bases on the RR 3 state. Additionally, we 
find two topologically-degenerate groundstate sectors on the 
infinite cylinder, respectively corresponding to the identity 
and the Eibonacci anyonic quasiparticle, which fully confirms 
the RR 3 state, without input of any features (such as shift) 
taken from the model wavefunction, that might have biased 
the calculation. The current work opens up a number of direc¬ 
tions deserving further exploration. Eor example, while the 
EQH ly = 12/5 state has been observed in experiment, there 
is no evidence of a EQH phase at ^ = 13/5 in the same sys¬ 
tems [5, 8 ]. So far it is not clear whether this absence is due to 
a broken particle-hole symmetry from Landau level mixing, 
or other asymmetry effects such as differences in the quantum 
wells [7]. Our numerical studies suggest that the outlook for 
the existence of such a state at 13/5 is promising, and some 
positive signs of this may have already been observed very re¬ 
cently [70]. Numerical studies may also further suggest how 
various other exotic EQH states in the second Landau level at 
different filling-factors may be stabilized. 

Note added .— After the completion of this work, we be¬ 
came aware of overlapping results in Refs. [73]. 
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In this supplemental material, we provide more details of the calculation and results which were not given in the main text. 
In Sec. I, we briefly summarize the Haldane pseudopotentials in disk and spherical geometries used for calculation in the main 
text. In Sec. II, we give a detailed analysis of edge excitations of the fermionic Read-Rezayi (RR) k = 3 state, based on the root 
configurations. In Sec. Ill, we introduce the effect of finite layer-width and show the evolution of entanglement spectrum (ES) 
with the change of the layer width. In Sec. IV, we extract the topological entanglement entropy (TEE) based on the dependence 
of entropy on the length of the orbital cut. In Sec. V, we show the ES of the = 12/5 state, which is a particle-hole-conjugate 
state of the z/ = 13/5 state. In Sec. VI, we introduce the numerical details of the infinite-size density-matrix renormalization 
group (DMRG) algorithm on cylinder geometry. 


I. Pseudopotentials 

Haldane[26] first pointed out that, any rotationally-invariant two-body interaction can be completely described by a set of 
“pseudopotential” V(m) with m > 0, if projected onto a single Landau level. The pseudopotential describes the energy of a pair 
of particles in a state of given relative angular momentum m. This formalism turned out to be useful not just for describing the 
details of the interaction, but also for understanding the microscopic conditions of the fractional quantum Hall (EQH) states in 
such systems. Here we focus on the pseudopotential formalism for two-body Coulomb interaction, in two specific geometries, 
the disk (or plane) and the sphere, respectively. 


1. Disk (plain) geometry 

In disk geometry, the Haldane pseudopotentials have been obtained many places [26]: 

CXD 

/ 2 

0 

Here |n, m) is a two electron state with relative (azimuthal) angular momentum m in the n-th Landau level. {x) is a Laguerre 
polynomial, and the two-body interaction V(r) = J (iqV(g)e®‘* ''. 

Eor the ideal Coulomb interaction in two dimensions V{r) = l/r (in unit of ^), we have, for the first Landau level, 

v("=°)(m)= 

2 m! 

and for the second Landau level: 

y(n=i)f , ^ r(m -f 1/2) (m - 3/8)(m - 11/8) 

2m! (m — l/2)(m — 3/2) 


2. Sphere geometry 


In spherical geometry, we first define the total angular momentum L = Q + n, where Q is the strength of the magnetic 
monopole in the center of the sphere and n is Landau level index. The pseudopotential V”(/) is defined as the interaction energy 
of a pair of electrons as a function of their parr angular momentum 1. In this expression, the pseudopotential for particles in a 
single Landau level is evaluated from the matrix element 

V”(0 =< n,L-l\V{ri - r 2 )\n,L-,l > 


Lor the Coulomb potential on the sphere, we define the chord distance between two points on a sphere as 


V{ri - r 2 ) = V{\ri - r 2 \) 


1 

i?V2 — 2 cos 012 


1 

R 


'^Pn{cOS 6 > 12 ), 


where Pn{x) is a Legendre polynomials. We omit the detailed calculations of integrals here [? ? ], and just present the final 
result: 


v"(0 


.. 2L 

k=0 


I L L 
k L L 


L k L\ 
-Q 0 Q ) 


( 2 ) 
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where 


L k L 
-Q 0 Q 


is the Wigner 3j coefficient and 


I L L 
k L L 


is the Wigner 6j coefficient. 


II. Edge mode counting based on “root states” 

Here we analyze the counting rules of the edge spectrum in the different topological sectors of the fermionic Read-Rezayi 
(RR) /c = 3 state. For simplicity, our analysis below is based on the highest density “root configurations”[63, 64], which are also 
the only-surviving configurations in the “thin-torus” limit of the FQH effect[62]. These obey characteristic “fractional exclusion 
statistics” rules that constrain the number of particles allowed in a certain group of consecutive orbitals, and the “admissible 
configurations” that obey these rules are in one-to-one correspondence with the states of the of zero-energy eigenstates of the 
model Hamiltonian for which the RR states are the highest-density zero-energy states. For the the k = 3 v = 3/5 RR state, the 
rules are [64]: “not more than one particle in any orbital” and “not more than three particles in any five consecutive orbitals”. 

We first assume that the lowest edge-mode with AL = 0 relates to the quantum Hall system with an open right edge. 
For example, for root configuration “... 11100111|000000...” on the cylinder, “|” separates the cylinder into left and right 
subsystems. This is the exclusion-statistics analog of a “filled Dirac sea, with the particles moved as far to the left as possible, 
consistent with the exclusion rules, (this is essentially the the defining relation of a Virasoro-primary state). The edge mode 
excitations can be obtained by the rightwards rearrangements of the particles at the edge, increasing the momentum above its 
minimum value. The excited configurations must still obey the exclusion rules, and this simple rule gives the multiplicities (or 
“characters”) of the spectrum as a function of “momentum” (Virasoro level) relative to the “primary” “Dirac sea” state. This 
method gives a simple exclusion-statistics-based method for obtaining the multiplicities that agrees with the very different and 
more opaque methods of conformal field theory (CFT) based on construction of “Verma modules”. 

There are two different topological sectors for the parafermion RR 3 state. One has the root configuration “... 0111001110 ...” 
and the other has “... 1010110101...”. All possible edge excitations of “... 0111001110 ...” at AL < 3 are listed in Tables 
I and II, which relate to the partitions with 3n and 3n — 1 electrons, respectively. The Table III and IV show the results for 
“... 1010110101 ...”. 


TABLE I: In this table, we analyze the counting rule of the edge excitations in the “... 11100111001110|00 ...” sector, which has multiplic¬ 
ities 1,1, 3, 6,... at AL = 0,1, 2, 3,.... 

0 

II 

< 

AL = 1 

AL = 2 

AL = 3 

II 

<1 

11100111001110|0000 

11100111001101|0000 

11100111001100|1000 

11100111001011|0000 

11100110101101|0000 

11100111001100|0100 

11100111001010|1000 

11100111000111|0000 

11100110101011|0000 

11100110101100|1000 

11010110101101|0000 

11100111001100|0010 

11100111001010|0100 

11100110101100|0100 

11100111001001|1000 

11100111000110|1000 

11100110101010|1000 

11100110100111|0000 

11100101101011|0000 

11010110101011|0000 

11100110011100|1000 

11010110101100|1000 

1101011010110101101|0000 
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TABLE II: In this table, we analyze the counting rule of the edge excitations in the ... 111001110011|00 ... sector, which has multiplicities 
l,2,5,9,...atAL = 0,l,2,3, .... 


o 

II 

<1 

AL = 1 

AL = 2 

AL = 3 

111001110011|0000 

111001110010|1000 

111001110010|0100 

111001110010|0010 


111001101011|0000 

111001110001|1000 

111001101010|0000 

111001100111|0000 

110101101011|0000 

111001110001|0100 

111001101010|0100 

111001101001|1000 

111001100110|1000 

111001011010|1000 

110101101010|1000 

110101100111|0000 

110100110101101011|0000 


TABLE III: In this table, we analyze the counting rule of the edge excitations in the ... 101011010110101| ... sector, which has multiplicities 
1,3,6,13, ...atAL = 0,1,2, 3,.... 


o 

II 

>-1 

<1 

AL = 1 

AL = 2 

AL = 3 

101011010110101|0000 

101011010110100|1000 

101011010110100|0100 

101011010110100|0010 


101011010110011|0000 

101011010110010|1000 

101011010110010|0100 


101011010101101|0000 

101011010101100|1000 

101011010101011|0000 

101011001110011|0000 

101010110101101|0000 

101011010101100|0100 

101011010110001|1000 

101011010101010|1000 

101011001110010|1000 

101011010011100|1000 

101010110101100|1000 

101011010100111|0000 

101011001101011|0000 

100111001110011|0000 

101010110101011|0000 

011010110101101|0000 


TABLE IV: In this table, we analyze the counting rule of the edge excitations in the ... 1010110101100|00 ... sector, which has multiplicities 
1,2,5,10,... at AT = 0,1,2,.... 


o 

II 

>-1 

<1 

AL = 1 

(M 

II 

>-1 

<1 

AL = 3 

1010110100|0000 

1010110010|0000 

1010110001|0000 

1010110000|1000 


1010101100|0000 

1010101010|0000 

1001110010|0000 

1010011100|0000 

0110101100|0000 

1010101001|0000 

1001110001|0000 

1010100110|0000 

1010011010|0000 

1001101010|0000 

0110101010|0000 

0110011100|0000 

011010110101100|0000 

100111001110010|0000 


III. Results for finite-layer thickness 

In a two dimensional system, the ideal Coulomb interaction between electrons has V{q) = 1/q. The finite thickness in the 
normal direction of an experimental quantum Hall system modifies the short-distance part of the ideal 2D interaction, yielding 
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an effective “softer” electron-electron interaction. Here we include this non-zero thickness effect in the Coulomb interaction 
through the standard Fang-Howard model [? ]. The Fang-Howard model can faithfully describe two dimensional heterostructure 
in FQH experiments, which assumed that the charge distribution normal to the x-y plane takes the form of variational wave 
function 77 ( 2 ) = exp(— 2 / 6 ), where b is the parameter giving the effective width of wavefunction in 2 -direction. The 

effective electron-electron interaction is then written as [? ] 


, 18 + 9qb + 3q^b^ 

= g 8 ( 1 + ,!,)» ■ 


(3) 
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FIG. 5; Oribital ES of the groundstate for partition with 3n (top) electrons and 3n — 1 electrons (bottom), for different layer thickness parameter 
b = 0.0, b = 1.0 and b = 2.0. The calculation is based on the pseudopotential obtained from plain geometry. The system size is Ne = 24 and 
Ns = 37. The first four degeneracy pattern are labeled by numbers. 


To study the effect of the finite-layer thickness, we use the pseudopotentials obtained from the infinite planar geometry (Eq. 
1) where the finite-layer thickness effect can be more conveniently obtained, although we study the quantum Hall systems on the 
sphere geometry. As noted before [65], the pseudopotentials in the spherical geometry approach those in the planar geometry if 
the spherical radius is taken to infinity in the thermodynamic limit. To support the above statement, we present the entanglement 
spectrum (ES) at 6 = 0, as shown in Eig. 5. The countings in the first several momentum sectors match the prediction from RR 
fc = 3 state. The similar picture obtained from the pseudopotential from planar geometry also indicates that the ES is robust 
and insensitive to the details of the pseudopotential form. We further show the result of nonzero layer thickness, in Eig. 5. The 
counting for the first four momentum sectors {AL < 3) are 1,1, 3, 6 for partition with 3n electrons in the subsystem. Examining 
the ES for different layer thickness, it is found that the ES deviates from the expected counting in AL = 3 around b Ri 2.0. 


IV. Topological Entanglement Entropy 

Eor a two-dimensional gapped topologically-ordered state, the entanglement entropy Sa{Ia) for a subsystem A with a finite 
boundary length I a is given by Sa{Ia) = c^Ia — 7 , where the TEE 7 is related to the total quantum dimension by 7 = lnT> 
[34, 35]. Since 25 contains the information about the quasiparticle content, the TEE can determine whether a given topological 
phase belongs to the universality class of a given topological field theory. 

Eig. 6 (a) shows numerically-calculated orbital-cut entanglement entropy S{Ia, Ng) as a function of the number of orbital (Ia) 
in the northern hemisphere for different system sizes (Ng). The initially-increasing parts of S{Ia, Ng) reflect the physics of the 
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FIG. 6: (a). Entanglement entropies with orbital partitioning S{lA,Ne), for various system sizes Ne, where Ia < ^Ne is the number of 
orbitals abive the cut in the northern hemisphere, (b). Scaling of entanglement entropies to S{Ia, Ne = 36) = —7 + ay/lA based on the solid 
circles (O the sphere, The length of the cut is proportional to ^/Ta [? ]). The dashed blue line is the theoretical RRi ;=3 TEE value 7 « 1.447. 
The open circles were discarded in the extrapolation because they represent very small subsystems (Ia < 5) and violate the area law [? ? ] 
and the finite-size saturation effect {Ia > 16). 


macroscopic state, while the downward curvature is a finite-size effect[? ? ? ]. With the help of the DMRG, we can obtain 
reliable entropies for Ia < 16, because the entropy for a given Ia < 16 is nearly saturated as Ng increases from 30 to 36, as 
shown in Fig. 6 (a). In Fig. 6 (b), we extract the TEE (red line), based on the raw data from Ng = 36. The obtained TEE is 
7 = 1.491 ±0.091 (If we perform the extrapolation based on different system sizes, similar results are obtained; for example, for 
Ne = 33, we get 7 m 1.519 ± 0.081). This is consistent with the theoretical value lnT> = In •\/5(l ± (j)) ~ 1.447 for RR k = 3 
state, where each non-Abelian Eibonacci anyon quasiparticle sector contributes quantum dimension dp = (p = (y/S ± l)/2 (cj) 
denotes the Golden Ratio). The appearance of c?f = ^ is a signal of the emergence of Eibonacci anyon quasiparticles, indicating 
two Eibonacci quasiparticles may fuse into the identity or into one Eibonacci quasiparticle. This exotic property makes Eibonacci 
quasiparticles capable of universal quantum computation, where all the quantum gates can be operated and measured by braiding 
Eibonacci anyons [17]. 


V. Entanglement spectrum at = 12/5 

Here we show the entanglement spectrum for the groundstate at f = 2 ± 2/5, which is the particle-hole (PH) conjugate state 
of 13/5 state. In spherical geometry, the highest density “root configuration” for the PH-conjugate RR fc = 3 state has a pattern 
of “0001100011000... 11000”. Consequently, there are also two distinct ways of partitioning, with 2n — 1 (labeled as P[l|l]) 
and 2n (P[0|0]) electrons in the subsystem. In Pig. 7, we show the ES of two partitions for Ne = 20. Por P[l|l] the leading ES 
displays the sequence of multiplicities 1,1,3, 6,12 in the first five momentum sectors AL = 0,1, 2,3,4. Por P[0|0], ES shows 
the multiplicity pattern 1,2, 5, 9 in the AL = 0,1, 2,3 momentum sectors. The ES picture found is exactly the same as that 
of the p = 13/5 state, taking into account PH conjugation so that P[0|0] (P[l|l]) for v = 12/5 relates to P[l|l] (P[0|0]) for 
p = 13/5. The low-lying ES structure is that predicted by the SU ( 2)3 CPT for the RR fc = 3 state, and provides a fingerprint of 
the topological order for the groundstate at p = 12/5. 


VI. Infinite DMRG on cylinder geometry 

The main results in this paper have been obtained with finite DMRG calculations in spherical geometry. Spherical geometry 
simulations are efficient for calculating the groundstate energy and related neutral gap. The corresponding ES can also be 
obtained. Nevertheless, one drawback of spherical geometry is that one needs to select a “shift” value S in the calculation. 
Usually, this value is determined by some empirical knowledge of the model wavefunction or one needs to compare results 
using different shifts. Another disadvantage of spherical geometry is that the sphere has genus zero so that it is not suitable for 
discussing the topological degeneracy for topological ordered state. 

An alternative strategy is to treat the cylinder geometry using ithe infinite DMRG algorithm [44, 69? ]. Here, we briefly 






11 


Ne=20, Ns=52, v=2+2/5 


10 


- 1 - 1--t- — J_ “ T -1-h- 1 0 n--1- _1_ _!_ -I- 



9 


8 


c 


5 


_3 


2 


6 


= 1 1 


6 


P[1|1] 


P[0|0] 


4 


I ■ I ■ I ■ I 

108 112 116 120 


4 


136 138 140 142 144 146 


FIG. 7: Orbital ES of the GS for partition and PIOlO] for FQH v — 12/5 state. The system size is Ne. = 20 and Ns = 52. The first 

four multplicity-patterns are labeled by numbers. 

introduce our implementation of infinite DMRG. In the infinite-size DMRG, we first start from a small system size. Then we 
insert several orbitals (for v = 13/5, we add ten orbitals each time) in the center, and optimize the energy by sweeping over 
the inserted orbitals. After the optimization, we absorb the new orbitals into the original existing system (here 5 orbitals are 
added to the left system and 5 to the right environment) and get the new boundary Hamiltonian. We repeat these insertion, 
optimizing and absorption procedures until both energy and entropy convergence are achieved by keeping a large number of 
states (M). There are some advantages of infinite DMRG over the finite DMRG. Compared with the finite DMRG simulation, 
infinite DMRG grows the system by several orbitals at each iteration and only sweeps the inserted part, thus the computational 
cost is significantly reduced. In the infinite DMRG algorithm, we do not need to set a “shift” S in the calculation. One can 
also access different topological sectors by randomizing the initial DMRG process. Details of the realization of infinite DMRG 
on cylinder geometry and the related benchmark for model Hamiltonian will be given elsewhere. Here we present some details 
related to the results shown in the main text. 

Working on the cylinder geometry, we choose the Landau gauge A = (0, Bx), which conserves the y momentum around the 
cylinder. The single electron orbitals in N-th Landau level are: 



(4) 


where Xj = j. 


fere Xj = = 1, 2,..., is the center in x axis and I is the magnetic length. Hn{x) is the Hermite polynomial. 

Introducing the destruction (creation) operator j) for 'tpNj, the Coulomb interaction can be written as 



(5) 


where the Coulomb matrix elements are 



( 6 ) 


If we only consider the second Landau level (setting Ni = 1), we can rewrite the Hamiltonian as 

V{m,n)a'lai+„ai+m+naj^rn+2n + 



(7) 


I n>0,m>0 


where V{m, n) is the matrix element derived from the modified Coulomb interaction. In this work, when studying on cylinder 
geometry, we choose the form of the modified Coulomb interaction as [69] 



V{ri,r2) 


( 8 ) 
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This is the form suitable for DMRG calculations on the cylinder. In the implementation, we kept all Coulomb interaction 
terms \V{m,n)\ > 10“® within the truncated range n < 4^,m < ^Lyl2. We have checked that the physical quantities 
remain qualitatively unchanged when the truncation range is varied. Here we would like to point out that one of the important 
advantages of the infinite DMRG is that it can easily deal with this type of the Coulomb interaction in the cylinder geometry. 
In traditional finite-size DMRG in cylinder geometry, an additional one-body potential U (x) is needed to avoid the electrons 
becoming trapped at the two ends of the finite cylinder, when studying the systems with Coulomb interaction between electrons. 
The infinite DMRG naturally overcomes this issue since it can access the actual results near the center by sweeping and edge 
effect should be suppressed when the length of the cylinder grows long enough to reach the fixed point for the state on the infinite 
cylinder. 

To access the topologically different groundstates on cylinder using infinite DMRG, we repeatedly start the infinite DMRG 
simulation for a given system size Ly for several times. In all calculations, we do not presume any empirical knowledge from 
model wavefunction. To be more explicit, we do not set a seed-state or an orbital configuration according to the root configuration 
in the initial DMRG process. We find that, in the relatively larger systems {Ly G [21, 24]), the system will automatically select 
one of the two groundstates jT'i) and j'l'^) with almost equal probability. In the smaller system size (Ly G [18, 20]), the system 
has larger probability to fall into jT'^) than |4'i). In all cases, we have double checked both of the groundstates are stable 
and robust with changing the parameters in DMRG calculation. Once one groundstate has been developed, the groundstate is 
robust against increasing keep states or increasing the cylinder length in DMRG. For example, in Fig. 8, we show the entropy 
evolution of one groundstate in jT*!) (red dots) and the other one in jT'^) (black dots). The datas come from two independent 
infinite DMRG simulations. At the point marked by M = 8000, we change the keep states from m = 7000 to to = 8000 
and keep to = 8000 for all steps after. The entropy of the two groundstates slightly increases with the increase of the keep 
states and the system length. There is no sign of tunnelling between the two groundstates in our infinite DMRG calculations. 
Physically, the tunnelling between two topological groundstates is forbidden since each topological groundstate hosts a well- 
defined anyonic flux line and the changing the global anyonic flux is energetically expensive. Furthermore, to check the two-fold 
groundstate degeneracy are complete, we try to start several infinite DMRG simulations with different random initializations. 
It is found that all simulations will randomly fall into or Thus we confirm the two-fold groundstates are indeed 

complete. Within randomizing the initial processes, although the two groundstates can be distinguished by characteristic orbital 
entanglement spectrum (as shown in main text), we And a small energy fluctuation and entropy fluctuation for each groundstate 
when we keep the same order of states in DMRG calculation. The entropy fluctuation is much smaller than the entropy difference 
AS Ri In^ Ri 0.47 between two topological sectors. We demonstrate the energy and entropy fluctuations as the error bars in 
Fig. 4(c-d) in the main text. The relatively larger uncertainty of entropy in system size Ly < 20 may result from; finite-size 
effect or that the groundstate at Coulomb point is very close to the phase boundary to a charge-density-wave (strip) state [25]. 
To clarify the latter possibility, we measure the mean orbital occupation number (uk) in the middle part of infinite cylinder. All 
of the groundstates have uniform occupation number with tiny fluctuation (Ank < 10“^). 
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FIG. 8: (Left) The entropy evolution of one groundstate realizing j'l',^) (black dots) and ['Ll) (red dots). These two groundstates come from 
two independent infinite DMRG simulations. We omit the infinite DMRG steps before we reach a nearly converged groundstate with keeping 
M = 7000. We change the keep state from M = 7000 to M = 8000 at the point as arrow marked. Each infinite DMRG step means inserting 
ten orbitals in the middle of the cylinder and sweeping until entropy converged. Here the system size is Ly = 24Zs. (Right) The mean orbital 
occupation number for two topological sectors on different system sizes. The dotted blue line shows v = 3/5. 




























